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NULL DEEPENING FOR AN ADAPTIVE ANTENNA BASED 
COMMUNICATION STATION 

FIELD OF INVENTION 

This invention relates to the field of wireless communication systems, and more 
5 specifically, to a method and apparatus for decreasing the transmitted power aimed at an 
interferer during transmission or the sensitivity to signals transmitted from the interferer 
during reception in order for a communications station to communicate to or from one or 
more remote subscriber units, the communication station having an adaptive antenna 
array and adaptive smart antenna processing. 

10 BACKGROUND TO THE INVENTION 

Adaptive smart antenna processing may be used in a communication station (e.g., 
a base station) equipped with multiple antennas to either reject interference when 
communicating from a subscriber unit to the communication station (i.e., on the uplink) 
or to deliver power in a spatially or spatio-temporally selective manner when 

15 communicating from the communication station to a subscriber unit (i.e., on the 

downlink). With smart antenna communication systems that use linear spatial processing 
for the adaptive smart antenna processing, during uplink communications, one applies 
amplitude and phase adjustments, typically but not necessarily in baseband to each of the 
signals received at the antenna array elements to select {i.e., preferentially receive) the 

20 signals of interest while minimizing any signals or noise not of interest — that is, the 
interference. Such baseband amplitude and phase adjustment can be described by a 
complex valued weight, the receive weight, and the receive weights for all elements of 
the array can be described by a complex valued vector, the receive weight vector. 
Similarly, the downlink signal is processed by adjusting the amplitude and phase of the 

25 baseband signals that are transmitted by each of the antennas of the antenna array. Such 
amplitude and phase control can be described by a complex valued weight, the transmit 
weight, and the weights for all elements of the array by a complex valued vector, the 
transmit weight vector. In some systems, the receive (and/or transmit) weights include 
temporal processing, and in such cases, the receive (and/or transmit) weights may be 

30 functions of frequency and applied in the frequency domain or, equivalently, functions of 
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time applied as convolution kernels. Alternatively, each convolution kernel, if for 
sampled signals, may itself be described by a set of complex numbers, so that the vector 
of convolution kernels may be re- written as a complex values weight vector, which, for 
the case of there being M antennas and each kernel having K entries, would be a vector of 
KM entries. 

Many methods are known for performing interference rejection or selective power 
delivery. Examples include least-squares beamforming and zero-forcing beamforming. 
Selective power delivery must balance competing goals. In general, the power delivered 
to one remote user cannot be simultaneously maximized while the power delivered to 
another remote user is minimized. More generally, if several remote users require power 
minimization (i.e., nulling), the relative power delivered to each must be traded off. This 
tradeoff can be based on a number of factors. For example, for a given interferer (which 
might be a co-channel user), a deeper null (i.e., decreased transmitted power aimed at an 
interferer during transmission or decreased sensitivity to signals transmitted from the 
interferer during reception) may be required for a remote user co-participating in the 
spatial channel established at the base station than for a remote user communicating with 
a different base station. 

A disadvantage of known methods, such as least-squares and zero-forcing, is their 
inability to flexibly perform this tradeoff. Zero-forcing methods attempt to direct perfect 
nulls towards all interferers, regardless of the power delivered to the desired remote user. 
Least-squares methods minimize a cost function which, when the uplink strategies are 
used to determine downlink strategies, on the downlink directs strongest nulls to remote 
users who were received most strongly on the uplink. Neither zero-forcing nor 
least-squares behavior may be appropriate for some systems. 

A particular case where least-squares nulling behavior is undesirable occurs when 
subscriber units are subject to the "near/far problem." In the near/far problem, one 
subscriber unit, say one denoted SU1, is far from the base station while a second 
subscriber unit, say one denoted SU2, is close. When the downlink weight vector, 
denoted wl, is computed for SU1, signals from SU2 are seen as interference and a null is 
formed towards SU2. The depth of the null in the direction of SU2 is limited by noise 
and other factors, including the possibility that the transmitted powers of subscriber units 
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are adjusted such that equal power levels are received by the base station. If high power 
is used to reach SU1 during base station transmit, then excessive signal levels leaking 
"through" in the imperfect null can disturb SU2. A deeper null towards SU2 is therefore 
desirable. 

5 A disadvantage of zero-forcing methods is the need for the base station to have 

complete knowledge of all remote user and interferer spatial signatures (or spatio- 
temporal signatures). The receive spatial signature and the receive spatio-temporal 
signature characterizes how the base station array receives signals from a particular 
subscriber unit in the absence of any interference or other subscriber units. The transmit 

10 spatial signature and the transmit spatio-temporal signature of a particular remote user 
characterizes how the remote user receives signals from the base station in the absence of 
any interference. See US patent 5,592,490 entitled SPECTRALLY EFFICIENT HIGH 
CAPACITY WIRELESS COMMUNICATION SYSTEMS, to Barratt et al, assigned to the 
assignee of the present invention and incorporated herein by reference, and US patent 

15 5,828,658 entitled SPECTRALLY EFFICIENT HIGH CAPACITY WIRELESS 

COMMUNICATION SYSTEMS WITH SPATIO-TEMPORAL PROCESSING, to Ottersten 
et al., also assigned to the assignee of the present invention and incorporated herein by 
reference, for spatial processing and spatio-temporal processing methods that use spatial 
and spatio-temporal signatures. Note that because a signature may be a spatial signature 

20 or a spatio-temporal signature, depending on whether the smart antenna processing is 
spatial or spatio-temporal, the term signature will be used herein, and whether the 
signature is spatial or spatio-temporal will depend on whether the processing is spatial or 
spatio-temporal, and whether the signature is a transmit or a receive signature will 
depend on the context, and which signature will be clear to those of ordinary skill in the 

25 art from the context. 

Determining all spatial signatures of all remote users and all interferers may often 
be impossible to carry out accurately if any of the interferers are weak or have a signal 
structure that is a priori unknown, or because of computational power limitations. 

In the uplink direction, increasing the depth of a null (i.e, decreasing sensitivity to 
30 signals received from a particular interferer) is desirable when the uplink strategy, for 
example the uplink weight vector, computed from a previous burst is used in predictive 
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mode, e.g., for new data. When uplink remote users are executing power control, for 

example to control received power at the base station or to reduce transmit power during 

o 

periods of voice or data inactivity, their uplink power may vary widely from burst to 
burst. Similar effects occur in fading environments. Thus, if a least-squares approach is 
5 used, the null depth obtained from a strategy computed for the previous burst may be 
inappropriate for the current burst. 

Thus there is a need in the art for a flexible method for directing precise and deep 
nulls on the uplink or downlink direction. There also is a need in the art for a method and 
apparatus for estimating one or more signatures in the direction of one or more interferers 
10 to use, for example, for such null deepening. There also is a need for a method for using 
a signature estimate of an interferer for directing precise and deep nulls in the direction of 
the interferer. There also is a need in the art for a method for directing precise and deep 
nulls in the direction of one or more interferers while substantially maintaining the other 
nulling and gain patterns of a provided adaptive smart antenna processing strategy. 

15 SUMMARY 

An advantage of the present invention is that it provides a flexible method for 
directing precise and deep nulls (i.e., decreasing the transmitted power aimed at an 
interferer during transmission or decreasing the sensitivity to signals transmitted from the 
interferer during reception) on the uplink or downlink direction. 

20 Another advantage of the present invention is that it provides a method and 

apparatus for estimating one or more signatures in the direction of one or more 
interferers. 

Another advantage of the present invention is that it provides a method for using 
a signature estimate of an interferer for directing precise and deep nulls in the direction of 
25 the interferer. 

Another advantage of the present invention is that it provides a method for 
directing precise and deep nulls in the direction of one or more interferers while 
substantially maintaining the other nulling and gain patterns of a provided adaptive smart 
antenna processing strategy. 
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Another advantage of the present invention is that it provides a null deepening 
method that essentially requires only an estimate signature of any interferer to deepen the 
null to that interferer, without needing a full set of signatures, and in particular, without 
needing the signature of the desired remote user. 

One asfpect of the invention is a method for determining, in a communication 
station using multiple antennas, improved uplink or downlink processing strategies, for 
example in tBe form of uplink or downlink weights for linear mart antenna processing, 
for which one or more nulls have a controlled depth. The method can be applied as a 
modification to a variety of known techniques for uplink and downlink strategy 
computation. The method needs as side information only the signatures of those remote 
users to wftich controlled nulls are to be directed. Other nulls in the array pattern, for 
which no signature estimates may be available, are substantially preserved. Another 
aspect of/the invention includes estimating the required signatures. Another aspect of the 
invention modifies existing uplink and downing strategy computation methods by 
15 injecting into the received signal sequence a synthetic signal whose signature equals that 
of the estimated remote user direction, and whose power is proportional to the desired 
null djpth. 

To overcome the limitations on null depth imposed by having a limited number of 
array samples, in an improved version of the invention, the signatures in the directions to 
20 be nulled are estimated by combining measurements over several bursts. An aspect of the 
invention identifies a particular type of signature estimation technique that combines 
favorably with the synthetic signal-injection method. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention will be more fully understood from the description of the 
25 preferred embodiments of the invention, which, however, should not be taken to limit the 
invention to any specific embodiment but are for purposes of explanation and to provide a 
better understanding of the inventive apparatus and method. The preferred embodiments 
may be better understood with reference to the following figures: 

Fig. 1 is a functional block diagram of a multi-antenna transceiver system which 
30 includes elements suitable for implementing the method of the present invention; 



Fig. 2 is a more detailed block diagram of a transceiver that includes a signal 
processor capable of executing a set of instructions for implementing the method of the 
present invention; 

Fig. 3 is a block diagram showing a strategy computation method and apparatus 
having as input a set of input signal data and a reference signal for the set of input data; 

Fig. 4 is a block diagram illustrating a reference signal generator that may be used 
in any of the embodiments of the invention; 

Fig. 5 is a block diagram illustrating a strategy computation method and apparatus 
that uses estimates of one or more characteristic features of a set of signal data; 

Fig. 6 is a block diagram illustrating some of the signal processing operations 
applied to received data when implementing an embodiment of the inventive method 
using the strategy computation of Fig. 3; 

Fig. 7 is a block diagram illustrating some of the signal processing operations 
applied to received data when implementing an embodiment of the inventive method for 
deepening more than one null; 

Fig. 8 is a block diagram illustrating some of the signal processing operations 
applied to received data when implementing an embodiment of the inventive method 
using the strategy computation of Fig. 5 with the characteristic feature being the 
covariance; 

Fig. 9 is a block diagram illustrating some of the signal processing operations 
applied to received data when implementing yet another embodiment of the inventive 
method using a noise-plus-interference covariance estimator; 

Fig. 10 is a block diagram showing the signal injection based embodiment of the 
inventive method applied to a non-linear processing strategy; 

Fig. 11 is a block diagram showing generalized combining method with the signal 
injection embodiment of the inventive method applied to a processing strategy; 
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Fig. 12 is a block diagram showing using a signal injection embodiment of the 
inventive method applied to a processing strategy generator to generate parameters that 
are then used to process downlink data; 

Fig. 13 is a block diagram showing a covariance modification embodiment of the 
5 of the inventive method applied to a covariance based processing strategy generator 
which uses a reference signal, primary data, and a covariance estimator; 

Fig. 14 is a block diagram showing an alternate embodiment of the combining of 
the input signed covariance estimate and interferer signature covariance estimate; 

Fig. 15 is a block diagram showing a covariance modification embodiment of the 
10 of the inventive method applied to a covariance based processing strategy generator 
which uses a reference signal, primary data, and a noise-plus-interference covariance 
estimator; and 

Fig. 16 is a plot of the performance obtained by using aspects of the invention in 
a simulated environment. 

15 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Adaptive Smart Antenna Processing 

The invention is directed to a processing method for altering the transmit or 
receive weights used by a communication station to define a transmitted signal or to 
process a received signal in order to deepen or otherwise manipulate the depth of a null 

20 formed to mitigate the effects of one or more known interferers. The interferer may or 
may not be another remote user sharing the same communication channel with the same 
base station. The invention may be implemented in a communication station that includes 
a receiver, an array of antennas and means for adaptive smart antenna processing of 
received signals. The invention may also be implemented in a communication station that 

25 includes a transmitter, an array of antennas, and means for adaptive smart antenna 

processing of transmitted signals. In a preferred embodiment, the communication station 
includes a transceiver and the capability of implementing both uplink and downlink 
adaptive smart antenna processing. 
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When receiving a signal from a subscriber (remote) unit, the signals received by 
each of the antenna array elements are combined by the adaptive smart antenna 
processing elements to provide an estimate of a signal received from that subscriber unit. 
In the preferred embodiment, the smart antenna processing comprises linear spatial 

5 processing, wherein each of the complex-valued (i.e., including in-phase / and quadrature 
Q components) signals received from the antenna elements is weighted in amplitude and 
phase by a weighting factor and the weighted signals are then summed to provide the 
signal estimate. The adaptive smart antenna processing scheme (i.e., the strategy) can 
then be described by a set of complex valued weights, one for each of the antenna 

10 elements. These complex valued weights can be described as a single complex valued 
vector of M elements, where M is the number of antenna elements. Thus, in the linear 
case, the smart antenna processing is designed to determine a set of weights such that the 
sum of the products of the weights times the antenna element signals provides an estimate 
of the remote user's transmitted signal which satisfies some prescribed "estimation 

15 quality" measure. 

This representation of the adaptive smart antenna processing can be extended to 
include spatio-temporal processing, where the signal at each antenna element, rather than 
being weighted in amplitude and phase, is filtered by a complex valued filter, typically for 
purposes of time equalization. In such a method, each filter can be described by a 
20 complex-valued transfer function or convolving function. The adaptive smart antenna 
processing of all elements can then be described by a complex valued M- vector of M 
complex valued convolving functions. 

Several methods are known for determining the weighting vectors to be applied 
when processing received signals. These include methods that determine the directions of 

25 arrival of signals from subscriber units, and methods that use the spatial or spatio- 
temporal characteristics of subscriber units, for example, the spatial or spatio-temporal 
signatures. See for example United States patents 5,515,378 and 5,642,353, entitled 
"SPATIAL DIVISION MULTIPLE ACCESS WIRELESS COMMUNICATION 
SYSTEMS", to Roy, et aL, assigned to the assignee of the present invention and the 

30 contents of which are incorporated herein by reference, for methods that use directions of 
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arrival. See also the above-referenced United States patents 5,592,490 and 5,828,658 for 
methods that use spatial and spatio-temporal signatures. 

"Blind" methods determine the weights from the signals themselves, but without 
resorting to a priori knowledge such as training signals or silent periods, that is, without 
5 determining what weights can best estimate a known symbol sequence (or in the case of 
the period silence, the absence of a known sequence). Such blind methods typically use 
some known characteristic of the signal transmitted by the subscriber unit to determine 
the best receive weights to use by constraining the estimated signal to have this property, 
and hence are sometimes referred to as property restoral methods. 

10 Property restoral methods in turn can be classified into two groups. Simple 

property restoral methods restore one or more properties of the signal without completely 
reconstructing the modulated received signal, for example by demodulating and then re- 
modulating. More complex restoral methods typically rely upon reconstruction of the 
received signal. 

15 Property restoral methods determine a signal (a "reference signal") that is 

constrained to the required property and then determine a set of weights corresponding to 
the reference signal, such that if the reference signal was transmitted by a remote user, the 
signals at the antenna elements of the receiving array would be acceptably "close" to the 
signals actually received. One example of a simple restoral method is the constant 

20 modulus (CM) method, which is applicable to communication systems that use a 

modulation scheme having a constant modulus, including, for example phase modulation 
(PM), frequency modulation (FM), phase shift keying (PSK) and frequency shift keying 
(FSK). The CM method has also been shown to be applicable to non-CM signals. Other 
partial property restoral techniques include techniques that restore the spectral properties 

25 of the signal, such as the signal's spectral self-coherence. 

"Decision directed" (DD) methods construct a reference signal by making symbol 
decisions (e.g.,, demodulating) the received signal. Such decision directed methods use 
the fact that the modulation scheme of the transmitted subscriber unit signal is known, 
and then determine a signal (a "reference signal") that is constrained to have the 
30 characteristics of the required modulation scheme. In such a case, the reference signal 
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production process includes making symbol decisions. Weights are determined that 
produce a reference signal, that if transmitted by a remote user, would produce signals at 
the antenna elements of the array that are acceptably "close" to the signals actually 
received. See, for example, United States patent application serial no. 08/729,390, 
5 entitled "METHOD & APPARATUS FOR DECISION DIRECTED DEMODULATION 
USING ANTENNA ARRAYS & SPATIAL PROCESSING" to Barratt, et ai, and serial no. 
09/153,1 10, entitled "METHOD FOR REFERENCE SIGNAL GENERATION IN THE 
PRESENCE OF FREQUENCY OFFSETS IN A COMMUNICATION STATION WITH 
SPATIAL PROCESSING" to Petrus, et al., both of which are assigned to the assignee of 
= io the present invention and the contents of which are incorporated herein by reference, for 

y3 descriptions of systems that use decision directed weight determination methods. 

2 As previously mentioned, weight determining schemes also are known that use 

*S training data, that is, data whose symbols are known a priori. The training data (possibly 

|7j with a timing offset or frequency offset, or both applied) is then used as a reference signal 

15 to determine the smart antenna processing strategy {e.g., the weights). Therefore, 

i=rr 

Efj reference signal based methods include the case in which the reference signal includes 

M training data, the case in which the reference signal includes a signal constrained to have 

2 some property of the transmitted signal, and the case in which the reference signal 
includes constructing a signal based on making symbol decisions. 

20 Non-linear uplink and downlink processing strategies also are known. In the 

uplink direction, such methods typically include demodulation and act to determine an 
estimate of the symbols transmitted by a desired remote user from the set of signals 
received at the antenna elements of the communication station. One known example of 
such a processing scheme is based on a Viterbi algorithm using branch metrics. In this 

25 regard, it is noted that the present invention is not limited to linear spatial and spatio- 
temporal processing methods that include weight determining, but also is equally 
applicable to non-linear methods such as those based on Viterbi algorithms and branch 
metrics, which may not necessarily include determining weights. 

In theory, adaptive smart antenna processing permits more than one 
30 communication link to exist in a single "conventional" communication channel so long as 
the subscriber units that share the conventional channel can be spatially (or spatio- 
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temporally) resolved. A conventional channel includes a frequency channel in a frequency 
division multiple access (FDMA) system, a time slot in a time division multiple access 
(TDMA) system (which usually also includes FDMA, so the conventional channel is a 
time and frequency slot), and a code in a code division multiple access (CDMA) system. 
5 The conventional channel is then said to be divided into one or more "spatial" channels, 
and when more than one spatial channel exists per conventional channel, the multiplexing 
is called space division multiple access (SDMA). SDMA is used herein to include the 
possibility of adaptive smart antenna processing, both with one and with more than one 
spatial channel per conventional channel. 

10 Base Station Architecture 

The preferred embodiment of the inventive method and apparatus is implemented 
in a communication receiver, in particular, a Personal Handyphone System (PHS)-based 
antenna-array communication station (transceiver) such as that shown in Fig. 1, with M 
antenna elements in the antenna array. The PHS standard is described, for example, in the 

15 Association of Radio Industries and Businesses (ARIB, Japan) Preliminary Standard, 
Version 2, RCR STD-28 and variations are described in Technical Standards of the PHS 
Memorandum of Understanding Group (PHS MoU — see http://www.phsmou.or.ip ). The 
preferred embodiments of the present invention may be implemented in two versions of 
the communication station of Fig. 1, one aimed at low-mobility PHS system, with M=4, 

20 and another, aimed at a wireless local loop (WLL) system, with a variable number, with 
typically M=12. 

While systems having some elements similar to that shown in Fig. 1 may be prior 
art, a system such as that of Fig. 1 with elements 131 and 133 capable of implementing 
the inventive method is not prior art. Note that the present invention is in no way 
25 restricted to using the PHS air interface or to TDMA systems, but may be utilized as part 
of any communication receiver that includes adaptive smart antenna processing means, 
including CDMA systems using the IS-95 air interface and systems that use the common 
GSM air interface. 

In the system of Fig. 1, a transmit/receive ("TR") switch 107 is connected 
30 between an M-antenna array 103 and both transmit electronics 113 (including one or 
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more transmit signal processors 119 and M transmitters 120), and receive electronics 121 
(including M receivers 122 and one or more receive signal processors 123). Switch 107 is 
used to selectively connect one or more elements of antenna array 103 to the transmit 
electronics 113 when in the transmit mode and to receive electronics 121 when in the 
5 receive mode. Two possible implementations of switch 107 are as a frequency duplexer 
in a frequency division duplex (FDD) system, and as a time switch in a time division 
duplex (TDD) system. 

The PHS form of the preferred embodiment of the present invention uses TDD. 
The transmitters 120 and receivers 122 may be implemented using analog electronics, 

10 digital electronics, or a combination of the two. The preferred embodiment of receivers 
122 generate digitized signals that are fed to signal processor or processors 123. Signal 
processors 119 and 123 incorporate software and/or hardware for implementing the 
inventive method and may be static (always the same processing stages), dynamic 
(changing processing depending on desired directivity), or smart (changing processing 

15 depending on received signals). In the preferred embodiments of the invention, processors 
119 and 123 are adaptive. Signal processors 119 and 123 may be the same DSP device or 
DSP devices with different programming for the reception and transmission, or different 
DSP devices, or different devices for some functions, and the same for others. Elements 
131 and 133 are for implementing the method of the present invention for downlink and 

20 uplink processing, respectively, in this embodiment, and include programming 
instructions for implementing the processing methods. 

Note that while Fig. 1 shows a transceiver in which the same antenna elements are 
used for both reception and transmission, it should be clear that separate antennas for 
receiving and transmitting may also be used, and that antennas capable of only receiving 
25 or only transmitting or both receiving and transmitting may be used with adaptive smart 
antenna processing. 

The PHS system is an 8 slot time division multiple access (TDMA) system with 
true time division duplex (TDD). Thus, the 8 timeslots are divided into 4 transmit (TX) 
timeslots and 4 receive (RX) timeslots. This implies that for any particular channel, the 
30 receive frequency is the same as the transmit frequency. It also implies reciprocity, i.e., 
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the propagation path for both the downlink (from base station to users' remote terminals) 
and the uplink (from users' remote terminals to base station) is identical, assuming 
minimal motion of the subscriber unit between receive timeslots and transmit timeslots. 
The frequency band of the PHS system used in the preferred embodiment is 1895-1918.1 

5 MHz. Each of the 8 timeslots is 625 microseconds long. The PHS system includes a 
dedicated frequency and timeslot for a control channel on which call initialization takes 
place. Once a link is established, the call is handed to a service channel for regular 
communications. Communication occurs in any channel at the rate of 32 kbits per second 
(kbps), a rate termed the "full rate". Less than full rate communication is also possible, 

10 and the details of how to modify the embodiments described herein to incorporate less 
than full rate communication would be clear to those of ordinary skill in the art. 

In the PHS used in the preferred embodiment, a burst is defined as the finite 
duration RF signal that is transmitted or received over the air during a single timeslot. A 
group is defined as one set of 4 TX and 4 RX timeslots. A group always begins with the 
15 first TX timeslot, and its time duration is 8 x 0.625 = 5 msec. 

The PHS system uses (tc/4 differential quaternary (or quadrature) phase shift 
keying (7C/4 DQPSK) modulation for the baseband signal. The baud rate is 192 kbaud. 
There are thus 192,000 symbols per second. 

Fig. 2 is a more detailed block diagram of a transceiver that includes a signal 
20 processor capable of executing a set of instructions for implementing the method of the 
present invention. This is the version of the Fig. 1 system suitable for use in a low- 
mobility PHS system. In Fig. 2, a plurality of M antennas 103 are used, where M = 4. 
More or fewer antenna elements may be used. The outputs of the antennas are connected 
to a duplexer switch 107, which in this TDD system is a time switch. When receiving, the 
25 antenna outputs are connected via switch 107 to a receiver 205, and are mixed down in 
analog by RF receiver modules 205 from the carrier frequency (around 1 .9 GHz) to an 
intermediate frequency ("IF"). This signal is then digitized (sampled) by analog to digital 
converters ("ADCs") 209. The result is then down converted digitally by digital 
downconverter 213 to produce a four-times oversampled complex valued (in phase / and 
30 quadrature Q) sampled signal. Thus, elements 205, 209 and 213 correspond to elements 
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that might be found in receiver 122 of Fig. 1. For each of the M receive tirneslots, the M 
downcon verted outputs from the M antennas are fed to a digital signal processor (DSP) 
device 217 (hereinafter "timeslot processor") for further processing. In the preferred 
embodiment, commercial DSP devices are used as timeslot processors, one per receive 
5 timeslot per spatial channel. 

The timeslot processors 217 perform several functions, which may include the 
following: received signal power monitoring, frequency offset estimation/correction and 
timing offset estimation/correction, smart antenna processing (including determining 
receive weights for each antenna element to determine a signal from a particular remote 
10 user, in accordance with the present invention), and demodulation of the determined 
signal. The version of the uplink processing method of the invention as implemented in 
each timeslot processor 217 in the embodiment of Fig. 2 is shown as block 241. 

The output of the timeslot processor 217 is a demodulated data burst for each of 
the M receive tirneslots. This data is sent to host DSP processor 231 whose main function 

15 is to control all elements of the system and interface with the higher level processing (i.e., 
processing which deals with what signals are required for communications in the different 
control and service communication channels defined in the PHS communication 
protocol). In the preferred embodiment, host DSP 231 is a commercial DSP device. In 
one implementation of the present invention, timeslot processors 217 send the determined 

20 receive weights to host DSP 231. Note that if desired, the receive weights may also be 
determined by software specifically implemented in host DSP 231. 

RF controller 233 interfaces with the RF transmit elements, shown as block 245 
and also produces a number of timing signals that are used by both the transmit elements 
and the modem. RF controller 233 receives its timing parameters and other settings for 
25 each burst from host DSP 231. 

Transmit controller/modulator 237 receives transmit data from host DSP 231. 
Transmit controller 237 uses this data to produce analog IF outputs which are sent to RF 
transmitter (TX) modules 245. The specific operations performed by transmit 
controller/modulator 237 include: converting data bits into a complex valued (ti/4 
30 DQPSK) modulated signal; up-converting to an intermediate frequency (IF); weighting by 
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complex valued transmit weights obtained from host DSP 231; and, converting the 
signals to be transmitted using digital to analog converters ("DACs") to form analog 
transmit waveforms which are provided to transmit modules 245. 

The downlink processing method of the invention is implemented in the 
5 embodiment of Fig. 2 in host DSP 23 1 , and is shown as block 243. In alternate versions, 
the downlink processing method is implemented in the timeslot processors 217 and in 
another version, it is implemented in transmit controller/modulator 237. 

Transmit modules 245 upconvert the signals to the transmission frequency and 
amplify the signals. The amplified transmission signal outputs are coupled to the M 
10 antennas 103 via duplexer/time switch 107. 

In describing the inventive methods, the following notation is used. Given M 
antenna elements (M=4 in one implementation, and 12 in another embodiment), let z\(t), 
Z2(t), . . ., z/ait) be the complex valued responses (that is, with in-phase / and quadrature Q 
components) of the first, second, M'th antenna elements, respectively, after down- 

15 conversion, that is, in baseband, and after sampling (four-times oversampling in the 

preferred embodiment). In the above notation, but not necessarily required for the present 
invention, t is discrete. These M time-sampled quantities can be represented by a single 
M- vector z(f) with the fth row of z(f) being z£t). For each burst, a finite number of 
samples, say N, is collected, so that z\(t) 9 z 2 (0, .. ., z^t) can each be represented as a N- 

20 row vector and z(f) can be represented by a M by N matrix Z. In much of the detailed 

description presented hereinafter, such details of incorporating a finite number of samples 
are assumed known, and how to include these details would be clear to those of ordinary 
skill in the art. 

Assume signals are transmitted to the base station from N s remote users all 
25 operating on the same (conventional) channel. In particular, assume that one of these, a 
particular subscriber unit of interest, transmits a signal s(t) t Linear adaptive smart antenna 
processing, which is used in the preferred embodiment of the invention, includes taking a 
particular combination of the / values and the Q values of the received antenna element 
signals z\{t\ Zii*), . . ., z^t) in order to extract an estimate of the transmitted signal s(t). 
30 Such complex valued weights may be represented by the receive weight vector for this 
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particular subscriber unit, denoted by a complex valued weight vector w r with I th element 
w Ti . The estimate of the transmitted signal from the remote unit may then be represented 

as: 

M 
/=1 

5 where w' xi is the complex conjugate of w ri and w J 1 is the Hermitian transpose (that is, the 
transpose and complex conjugate) of receive weight vector w r Eq. 1 is called a copy 
signal operation, and the signal estimate s(t) thus obtained is called a copy signal. 

The spatial processing described by Eq. 1 may be re-written in vector form for the 
case of N samples of M- vector signals z(t) and AT samples of the transmitted signal s(t) 
10 being estimated. In such a case, let s be a (1 by AO row vector of the N samples of s(t). 
The copy signal operation of Eq. 1 may then be re- written as s= w r Z . 

In embodiments which include spatio-temporal processing, each element in the 
receive weight vector is a function of time, so that the weight vector may be denoted as 
w r (r), with ith element w ri (t). The estimate of the signal may then be expressed as: 
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M 

s(t) = y £w' Ti (t)*z i (t) (2) 



where the operator "*" represents the convolution operation. Spatio-temporal processing 
may combine time equalization with spatial processing, and is particularly useful for 
wideband signals. Forming the estimate of the signal using spatio-temporal processing 
may equivalently be carried out in the frequency (Fourier transform) domain. Denoting 
20 the frequency domain representations of s(t\ z t {t) 9 and w n -(f) by S(k), Z t {k), and W/(fc), 
respectively, where k is the discrete frequency value: 

M 

S(*) = 5>ri (*) z /(*)- (3) 

t=l 

With spatio temporal processing, the convolution operation of Equation (2) is 
usually finite and when performed on sampled data, equivalent to combining the spatial 
25 processing with time equalization using a time-domain equalizer with a finite number of 
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equalizer taps. That is, each of the w n (r) has a finite number of values of t and 
equivalently, in the frequency domain, each of the Wj(k) has a finite number of k values. 
If the length of the convolving functions w n (f) is K, then rather than determining a 
complex valued M- weight vector w r , one determines a complex valued M by K matrix W r 
whose columns are the K values of w r (f). 

Alternatively, a spatial weight determining method can be modified for spatio- 
temporal processing according to a weight matrix by re-expressing the problem in terms 
of matrices and vectors of different sizes. As throughout this description, let M be the 
number of antenna elements, and N be the number of samples. Let K be the number of 
time equalizer taps per antenna element. Each row vector of N samples of the (M by AO 
received signal matrix Z can be rewritten as K rows of shifted versions of the first row to 
produce a received signal matrix Z of size (MK by AO, which when pre-multiplied by the 
Hermitian transpose of a weight vector of size (MK by 1), produces an estimated received 
signal row vector of N samples. The spatio-temporal problem can thus be re-expressed as 
a weight vector determining problem. 

For example, for covariance based methods, the weight vector is a "long" weight 
vector of size (MK by 1), the covariance matrix R ZZ =ZZ H is a matrix of size 
(MK by MK), and the correlation of the antenna signals Z with some signal represented 
by a (1 by N) row vector s is r zs =Zs H , a long vector of size (MA' by 1). Rearranging terms 
in the "long" weight vector provides the required (M by AO weight matrix. 

A downlink (i.e., transmit) processing strategy using adaptive smart antenna 
processing includes transmitting a signal, denoted in the finite sampled case by a (1 by AO 
vector s, from the communication station to a particular remote user by forming a set of 
antenna signals (typically, but not necessarily in baseband). Linear smart antenna 
processing determines the antenna signals as: 

Z = w t s, 

where w t is the downlink (or transmit) weight vector. When SDMA is used to transmit to 
several remote users on the same (conventional) channel, the sum of w u * for different 
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signals s f - aimed at different remote users is formed, to be transmitted by the M antenna 
elements. 

Note that a downlink strategy, for example including determining the downlink 
weights w t , may be implemented by basing it on an uplink strategy, for example uplink 

5 weights, together with calibration data. In this situation, the calibration accounts for 

differences in the receive and transmit electronic paths for the different antenna elements. 
The downlink strategy may also be found from the uplink strategy by using the transmit 
spatial signatures of the remote users, or by other known methods. Again, the weight 
vector formulation may be used for both linear spatial processing and linear spatio- 

10 temporal processing. 

Consider a communications station communicating with N s remote users on any 
channel. On the uplink, for an M antenna system, the M received signals at the base 
station can be stacked into an (Mby 1) received signal vector z(0 which may be modeled 
as 



where t is time, which in the preferred embodiment is discrete, Sj(r), /=1, N s , denotes 
the signal (preferably in baseband) at time t transmitted by the ith of N s remote users, aj 
is the spatial (or spatio-temporal) signature (a complex valued M- vector) of the ith 
remote co-channel user, and v(f) denotes additive noise (which may include other 
20 interfering signals, viewed as noise for any remote user). Note that the signature of any 
user is the signal received at the M antenna elements in the absence of noise and the 
absence of any other interfering users when the user transmits a unit impulse signal (see 
above referenced above-referenced United States patents 5,592,490 and 5,828,658). In 
matrix form the above signal model can be represented as 
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Ns 

z(o = X a / s ;W + v(r) ' 



(4) 



25 



z(0 = As(r) + v(f). 



(5) 



where 



A - t a i a 2 • a Afr]> ^ s(0 = bl (0*2 (0 - - *Ns (OF • 



(6) 



with the superscript T denoting the matrix transpose. 
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Processing Strategy Computation Methods 

The preferred embodiment of the invention improves a method for computing an 
uplink or downlink processing strategy that uses as inputs the received antenna signal 
data and typically a reference signal, and that takes into account the interference 
environment present in the received antenna signal data for interference mitigation. The 
improvement is to deepen or otherwise modify the null depth the strategy provides for 
mitigating the effects of any one or more known interferers. As would be known to those 
in the art, interference mitigating strategy determining methods explicitly or implicitly 
use one or more characteristic features of the received antenna data. Depending on the 
known strategy computation method, the input signal data may be explicitly reduced to 
one or more particular characteristic feature that the known method explicitly uses for its 
computation. For example, for methods that use the spatial or spatio-temporal covariance 
matrix of the input, the data may be reduced to the spatial or spatio-temporal covariance 
matrix of the data. Other methods may be based on other properties and in such cases the 
input signal data may be reduced to the particular feature or property which the known 
strategy method utilizes. Yet other methods, while implicitly dependent on some property 
(e.g., spatio-temporal covariance) of the received signal, do not require explicit 
estimation of the property such as the covariance. 

Since some of the discussion below applies for both receive and transmit strategy, 
the "r" or "t" subscripts are omitted in such quantities as the weight vectors w. Such 
subscripts may be used explicitly to identify uplink or downlink processing, and their 
addition will be clear from the context to those of ordinary skill in the art. 

Let Z be the matrix of received antenna array signals, preferably but not 
necessarily in baseband. Let Sj be a (1 by Afj) reference signal vector of samples. A 
reference signal s ls may be a known training sequence, or, in decision directed methods, a 
signal constructed to have the same known modulation structure as the signal transmitted 
by the subscriber unit, or for property restoral methods, a signal that is constrained to 
have the required property. 

The invention may be applied to deepen the nulls resulting from uplink or 
downlink strategies determined with any strategy computation method. Thus, in all of the 
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embodiments of the invention to be described herein, a known strategy (e.g., weight) 
determining process is used that computes an uplink or a downlink strategy based on 
received signal inputs Z from the antenna array elements. Many such methods also use 
reference signal, denoted S] determined for the remote user signal of interest from the , 
received signal data Z. Such a strategy determining method is shown in Fig, 3, which is 
marked prior-art, but is only prior art without the modifications described herein. The 
reference signal is typically extracted from a set Z of received data as shown in Fig. 4. 
Fig. 5 is one case of the system of Fig. 3 wherein the method explicitly uses as an input at 
least one of the characteristic features of the set of received signal data Z, so in Fig. 5, the 
set of input data is further operated upon to obtain one or more characteristic features 
(e.g., the covariance, the principal components of the covariance matrix, a particular 
feature of the data, etc.) of the data, which is then explicitly used in the known strategy 
(e.g., weight) determining method. 

One example of a weight determining method that has the structure shown in 
Fig. 5 is the well known least squares (MSE) technique computes the uplink or downlink 
weights by solving the minimization problem: 

w = argmin w |w H Z - s x || = (R zz r l r zs = (ZZH)-l Zsf 1 , (7 A) 

where R ZZ =ZZ H is the spatial (or spatio-temporal) covariance matrix of the antenna 
signals, || || denotes the vector norm, and r zs = Zs|* is the cross correlation between the 

antenna signals and the reference signal. The method of Eq. (7 A) is called the minimum 
mean squared error (MMSE) method herein, and uses the covariance R zz of the input 
signal data as the characteristic feature for weight determining. Thus, calculation of the 
weights according to Eq. (7 A) requires having data corresponding to the signals received 
at the antenna elements and a reference signal. In practice, this typically entails 
identifying the reference signal data and extracting it from a data burst, forming the 
covariance matrix for the received signals in the data burst, forming the cross correlation 
term, and solving for the weights. 

In practice, for example in a mobile PHS, because of limited computational 
power, only a small number of the 960 samples making up a received signal burst are 
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used to determine R zz _1 . In addition, the samples are determined at the baud points rather 
than being oversampled. Using only a small number of samples to determine weights may 
cause what is called "overtraining" herein: the weights performing well on the data 
containing the reference signal but not on new data, i.e., while the weights extract the 

5 desired signal energy and reject interference, the weights may perform poorly on new 
data. An improvement to the least squares technique which can assist in this situation 
includes a process termed "diagonal loading" (see, for example, B.L. Carlson: 
"Covariance matrix estimation errors and diagonal loading in adaptive arrays," IEEE 
Transactions on Aerospace and Electronic Systems, vol. 24, no. 4, July 1988), by which a 

10 diagonal adjustment is added as follows: 



where y is a small adjustable factor used to improve the performance of the least squares 
solution by reducing sensitivity to statistical fluctuations in Z. 

The methods described by Eq. (7B) also is of the structure shown in Fig. 5, again 
15 with the characteristic feature being the spatial (or spatio-temporal) covariance 
R zz = ZZ H . Other uplink processing methods are known that use the noise-plus- 
interference covariance matrix, denoted R vv , rather than the noise-plus-interference-plus- 
signal covariance matrix R zz . Such methods include some noise plus interference 
covariance estimator which uses the received data (containing the signal plus interference 
20 plus noise) and the reference signal for the desired user to determine R vv . 

The invention includes two main aspects: (i) modifying the uplink or downlink 
strategy to realize an improved null, and (ii) accurately estimating the signature in the 
direction where a modified (e,g., deeper) null is desired followed by null deepening 
based on the estimate. 

25 Null depth is generally limited by two effects. The first effect is the "natural" null 

depth of the uplink or downlink strategy. For example, when the signatures of the remote 
user and interferer are highly correlated, when the interferer power is not too strong, and 
when the signal-to-noise ratio (SNR) is modest, an uplink strategy based on maximizing 
the signal-to-interference-plus-noise ratio (SINR) will not direct particularly deep nulls 



w = (ZZ H + yI)- 1 Zsf 1 , 



(7B) 
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towards the interferer. One aim of the present invention is to deepen the nulls formed 
under these conditions. 

While a preferred embodiment is described in detail for deepening a null based on 
known or estimated signatures towards which the nulls are to be deepened, various 
embodiments are possible. They include using: synthetic signal injection, synthetic signal 
injection only for the "noise plus interference" estimation block, low-rank update of the 
noise plus interference estimate, and strategy orthogonalization. In the latter case, the 
uplink weights for two remote users that share a spatial channel are combined to 
construct weights with improved nulling performance, while again substantially 
preserving the rest of the null and gain pattern. These embodiments are detailed 
hereinunder. 

The described methods for null deepening use a known or estimated signature. 
Known signature estimation methods such as simple (e.g., maximum likelihood) 
signature estimation of the interfering remote user may not produce a sufficiently 
accurate signature estimate for effective null deepening in the direction of the interfering 
remote user. This is because the estimate may suffer from contamination by the signal of 
the remote user to which delivered power is maximized. In the preferred embodiment, 
the method used includes joint signature estimation using the reference signals of both 
the remote user and any interferer(s). Geometric methods based, for example, on angle of 
arrival, may alternatively be used to estimate signatures which then are used for null 
deepening according to any of the signature-based null deepening embodiments 
described herein. 

Another effect that limits null depth is the accuracy to which signatures or 
covariance matrices are estimated, which in turn is limited by the number of samples, the 
SNR of the remote user whose signature is to be estimated, and the power of other 
remote co-channel users. An additional aspect of the present invention is improving the 
signature estimate by combining data gathered over several bursts with reference signals 
of multiple remote users computed over the same bursts. This aspect of the invention is 
motivated by the recognition that the signature estimate of the interfering remote user 
must be accurate to a degree. To improve the null depth of a least-squares solution over a 
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single burst, for example, the signature estimate typically may need be based on estimates 
derived from several bursts. 

Null Deepening 

Using the invention can deepen null towards one remote user or more than one 
5 users based on knowledge of the signatures of the one or more users. The signatures may 
be known or estimated, and one aspect of the invention is a method for estimating the 
signatures applicable to null deepening. For best performance, the signature estimate 
should be substantially free of contamination from signals to which beams are directed. 
The invention uses the signature (known or estimated) of one or more known interferers 
10 to modify a known strategy or known strategy determining method to produce a modified 
strategy with improved null deepening. The knowledge of the signatures may be, for 
example, from historical signature records. 

For the description below, the desired remote user is denoted SU;, and the weight 
vector for communicating with SU; is denoted by w f . When communicating with SU /? all 
15 other users, denoted SUj where j is not i, are interferers, and the signature for such 
interferers are denoted ay, where j is not /. 

The Signal Injection Method 

A first embodiment of the invention applicable to one null-deepening based on 
the known or estimated signature of one interferer denoted SUj with signature (estimate) 
20 2ij is shown in Fig 6. A signal formed from the estimated or known interferer signal is 
scaled and added to the received signal data Z before determining the strategy, that is, 
modified received signals described by 

Z = Z +CLj&j& synt i l9 

with OLj a tunable scale factor for any interferer j. is any signal that has the same 
25 temporal structure as a typical communication waveform used in the system, or 
preferably, a temporal structure that would be identified by the particular strategy 
computation method as not that of the user. For example, for a constant modulus property 
restoral method, using a non-constant modulus structure for s^ nr ^ ensures that the 
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strategy computation method recognizes this as an interferer. In the preferred 
embodiment, however, comprises N random noise samples. In another 
embodiment, this is a constant signal, that is, the signal [1 1 ... 1]. The scale factor a,- for 
any interferer j can be set to a large fixed value to force a deep accurate null towards that 
5 undesired subscriber unit signature. The modified input signal Z is then applied to any 
strategy determining scheme, including the ones shown in Fig. 3 and Fig, 5. 

A generalization of this embodiment, called "signal injection" herein is illustrated 
in Fig. 7. Here, nulls are directed to several interferers denoted SUy where j is not i. 
Signals are formed from the estimated (or known) signatures of the undesired subscriber 
10 units, then are scaled and added to received signal data Z before determining the strategy, 
that is, modified received signals described by 

Z = Z+X a 7 a 7 S ^_7 > (8) 

with (Xj a tunable scale factor for any interferer j, are used in place of Z and applied to the 
strategy determining method, say that of Fig. 3 or that of Fig* 5. for processing SUi 

15 signals. The scale factor ot, for any interferer j can be set to a large fixed value to force a 
deep accurate null towards each undesired subscriber unit signature. Each of the S^^j 
is any signal that has the same temporal structure as a typical communication waveform 
used in the system and that would be recognized a non signal of interest, and Ssy, nt hj f° r 
any interferer j must be non-collinear with any s^^j* for any other interferer/. Again, 

20 in the preferred embodiment, each j is a set of random samples, with the sets for 
two distinct interferers being statistically independent. 

An improved embodiment makes ay a function of the relative power levels 
transmitted to SU f and yth subscriber units, or to any similar measure derived from 
received signal strengths, thus forming deep nulls in a controlled way only where needed. 
25 That is, the nulls are strong where strong nulls are needed and weak where such strong 
nulls are not needed. Thus forming weak nulls where allowable would reserve degrees of 
freedom for use elsewhere, for example to null additional noise sources. , 
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The parameter setting method is described in more detail later in this document. 

Fig. 8 shows an embodiment of the signal injection method for the case of null- 
deepening for one interferer for a co variance based strategy determining method. Fig! 9 
shows an embodiment of the signal injection method for the case of null-deepening for 
5 one interferer for an interference-plus-noise covariance based strategy determining 
method. 

Fig. 10 is a block diagram showing an embodiment of the inventive method 
applied to a non-linear strategy generator, where the strategy generator includes a 
demodulation stage. This may be a Viterbi-algorithm-based decoder which operates based 

10 on branch metrics, for example. The invention is shown supplying a modified noise-plus- 
interference covariance to a demodulator, with the modification based on the signal 
injection method. The resulting strategy is then applied to some new input data. The 
demodulator, in this case, may be a Viterbi-algorithm based demodulator in which 
covariance information is used to "whiten" {i.e., decorrelate) the input data, in this case 

15 the new input data. In an alternative, the reference data might also be used to estimate the 
channel, and provided in the form needed for the Viterbi-method in the demodulator. 

Note that while the signal injection methods preferably involve adding a fraction 
of the "synthetic" signal generated from the one or more known or estimated interferer 
signatures, other, non-additive, methods of combining also are within the scope of the 
20 present invention. This is shown in Fig. 11. For example, one method for combining 
includes performing a matrix factorization of the input data and the synthetic data into 
factors, and then combining the resulting factors to form a combined signal. The 
factorization may be a generalized singular value decomposition. 

Fig. 12 is a block diagram showing a general application of the inventive method 
25 as used to process data in the downlink data. As shown in Fig. 12, an uplink strategy 
generator {e.g., a weight determining method) is based on a reference signal (obtained 
from the primary data) and a set of data. In accordance with the present invention, the 
data input to the strategy generator is a combination of the input data and a signal 
generated from one or more known or estimated interferer signatures. The output of the 
30 strategy generator is a set of parameters which is used in the processing of the downlink 
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data. In the linear spatial processing case, the parameters may be a set of weights. Note 
that the parameters produced by the strategy generator may be combined with calibration 
data to produce the required downlink strategy parameters for processing of downlink 
data. 

5 Strategy Orthogonalization Method. 

According to this implementation, for the ith remote user, a strategy (e.g., a 
weight vector w , ) is determined that is orthogonal to the to estimated (or known) 
interferer signature(s) ay. Moreover, the strategy is determined from the strategy (e.g., a 
weight vector w £ ) for the desired user SU, and the strategy (e.g., a weight vector Wj) used 
10 to communicate with the interferer j. The strategies for the user and interferer, 

respectively, for example weight vectors w,- and w,, respectively, may be obtained using 
some known strategy determining method, for example that of Fig. 3 or Fig. 5 (including 
using Eq. (7 A) or Eq. (7B)). 

In the linear strategy case, assume a linear strategy w f (spatial or spatio-temporal) 
15 exists for receiving from or transmitting to a desired remote user, say the ith, and another 
linear strategy Wy for receiving from or transmitting to an interferer, say the 7th user, 
where j * i, where w,- and w ; have substantially the same null and gain patterns, except of 
course that w f - nulls the interferer and w y - nulls the desired remote user i because when 
communicating with the interferer, the remote user is an interferer. According to this 
20 embodiment of the invention, the improved strategy is 

W ( - = W/ + OLj Wy, (9) 

where the scalar tunable factor a, is chosen to make w ( - orthogonal to the interferer 
signature ay. That is, <Xy is chosen so that a^w,- = 0. This is obtained with 
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A generalization when one has a number of interferers to which nulls are to be 
deepened is, given a strategy w,- to a desired user t, and several weight vectors Wy, j * i, 
for interferers, then the strategy (i.e., weight vector) to use is 



where the OLj are chosen to ensure that the modified strategy w f is orthogonal to the 
interferer signatures apj ^ /. 

Another particular orthogonalization implementation that modifies the strategy w 
to be orthogonal to all known interferer signatures is,. 



where factor 7^ is chosen to make w, orthogonal to all the interferer signatures a^y * i. 
The factors jj may be computed using Gram-Schmidt orthogonalization, which in the 
case of a single interferer j leads to 



A disadvantage of the this forced signature orthogonalization method, however, is that it 
may disturb nulls which were formed towards "coherent" noise sources (e.g., from 
neighboring cells in a cellular system) or towards other incoherent interferers, and may 
have an effect on the "main lobe" towards the desired SU. Therefore, the preferred 
strategy orthogonalization method is that of Eq. (9) or its extension to multiple 
interferers. 

Covariance Modification 

A third embodiment is applicable to strategy determining methods that use an 
estimate of the spatial or spatio-temporal covariance (including the interference-plus- 
noise covariance) determined from the input signal (and possible a reference signal) to 
determine the strategy. The methods of Eqs. (7 A) and (7B) so use the covariance to 
determine a weight vector. One embodiment of the invention adds some information 



W; = W f + 
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from the interferer signatures to data used by the covariance estimator of the known 
weight determining method. For example, if the known weight computation method of 
Eq. (7B) is used, according to this embodiment of the invention, the weights to use are 
determined as: 



w = 



ZZ H +}l+£a 7 -a;a*J 



l 



ZsP, (10) 



where Z is the received data, and a,- again a tunable scale factor for any interferer j 9 set as 
described above for the signal injection method. 

Many modifications are possible within the scope of the invention. In general, any 
uplink (i.e., receive) processing method which takes as input a spatial covariance matrix 

R zz = ZZ H computed from input data Z can be modified by altering the spatial 

covariance matrix R zz to incorporate the effects of the interferers towards which the nulls 
are desired to be deepened, the modification being: 



ZZ H + 5>7 a 7 a ? 



(11) 



Note that in the preferred embodiment, the weight determining method to which the 
15 invention is applied takes as input a spatial covariance matrix R zz and a reference signal 

t_i 

cross correlation r zs = Zs £ * . 

Fig. 13 shows one version of this embodiment. The weight strategy computation 
method is one that uses covariances to determine spatial (or spatio-temporal) weights. 
These weights might be uplink weights, or may be downlink weights. The known scheme 

20 is modified using an aspect of the present invention to determine weights that include 
interference mitigation based on signature data on known interferers. The computation 
required includes estimating the known interferer spatial signatures when such spatial 
signature estimates are unavailable, and also estimating the covariance of those signature 
estimates and adding a tunable fraction of these interferer covariances to the covariance 

25 estimate determined from the received data. 
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Other covariance matrix modifications also are possible within the scope of the 
invention. In yet another alternative shown in Fig. 14, rather than the combining of the 
interferer data being by adding part of the interferer covariance matrix to the input signal 
covariance matrix, other methods of combining could be used. For example, the method 
5 for combining may includes performing a matrix factorization of the input data and of the 
signatures into factors, and then combining the resulting factors by one of several 
methods to form a combined covariance matrix. 

This embodiment can also be applied to processing methods that use the noise- 
plus-interference covariance matrix, denoted R vv , rather than the noise-plus-interference- 

10 plus-signal covariance matrix R zz . Such methods include some noise plus interference 

covariance estimator which uses the received data (containing the signal plus interference 
plus noise) and the reference signal for the desired user to determine R vv . An 
embodiment of the invention used in such a method is shown in Fig. 15, which shows 
one embodiment of the method of the invention applied to such R vv -based techniques, 

15 with R vv substituted for R zz . That is, for R vv obtained from the signal data, a 
modification 



f \ 
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R vv + X a 7 a 7 a ^ 

is made, where again the parameters 0Cj are tunable, and preferable determined from 
power control as described hereinbelow. 

20 Setting the Tunable Parameters 

Many strategies for selecting the parameters (Xy are possible within the scope of 
the invention, and some possible strategies are now described. 

Strategy 1 (applied, for example, to Signal Injection) 

Recalling that according to the signal injection embodiment, after obtaining 
25 precise signature estimates &j (or knowing the signatures) of the interferers, given a 
strategy computation method for determining the strategy for user i from the received 



* % 



30 

signal Z (and possible a reference signal for user 0, from the received signal Z is 
modified to 

Z = Z + S^Oty % s synthJ 
where s^^j is an artificially synthesized signal for the yth interferer. According to a 

5 first parameter determining strategy, the parameter ot, is a large constant chosen to be 
sufficiently large to ensure that the carrier to interference ratio (CIR) in the modified 
received signal Z is small. In this case, the base station is forced to place a deep null 
towards the jth subscriber unit, where j*i, for a strategy to the ith user. The carrier to 
interference ratio is easily estimated as the ratio of what is indicated by the received 

10 signal strength indicator (RSSI) typically available at a receiver to the interferer power, 
which is being controlled by a large value of the parameter. 

It is interesting to note that for an interference mitigating strategy computation 
method that uses a reference signal for interference mitigating, such as those of Eq. (7 A) 
and (7B), in the limit when (Xj — > <~, the transmit weight for the ith user is orthogonal to 

15 a,-. 

For the covariance modification method, note that the covariance of an input 
signal includes terms from each of the constituent signals in the input received signal Z. 
A practical covariance estimate may include cross terms, so approximately includes 
terms from each of the constituent signals in the input received. The covariance 
20 modification changes the interference constituents reflected in the covariance estimate. 
According to this strategy 1 when applied to the covariance modification method, the 
adjustable parameter for any interferer is selected to be a number sufficiently large to 
ensure that the carrier to interference ratio reflected by the modified covariance estimate 
is small. 

25 Strategy 2 (applied, for example, to Signal Injection) 

This strategy uses some knowledge of the transmitted power towards or from 
each subscriber unit. On the downlink, this information is usually readily available in a 
typical base station in the downlink power control unit. According to this parameter 
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tuning strategy, the parameters (Xj are chosen in such a way that the power distribution of 

the constituent signals in Z is inversely proportional to the distribution of the downlink 
power that the base station transmits towards the subscriber units. 

In order to find exact formulas for ay, one notes that the received signal Z may be 
5 modeled as 

Z = a/SfhS^^-Sy-h v, (13) 
where v is (non-co-channel) interference plus noise, a, and s z - are the signature of and 
transmitted signal from the remote user of interest, and the a y and sj are the signatures of 
and transmitted signals, respectively, of the interferers, where Recalling that 
10 according to the signal injection embodiment, after obtaining precise signature estimates 
&j (or knowing the signatures) of the interferers, the received signal Z is modified to 

Z = Z 4" ^tj^fiLj &j S^y n/ /j J 

where s^^j is an artificially synthesized signal for the yth interferer. Combining the 
model of Eq. (13) with the signal injection to give the following equation for the 
15 modified received signal: 

Z = afit+Lfri (afj +Oj afsynth j) + n - 
Assuming that the signals Sj and s^^j are uncorrected, the average powers of the 
constituent signals in Z due to the desired user / and the jth interferer, y*/, are given by 
la/lfo,- 2 and { ||a y - 2 aj 2 + lotyl 2 ||a 7 -|| 2 o 2 synthj}* respectively, where of denotes the 

20 average power of signal Sj and o 2 sy n thj denotes the average power of the synthesized 

signal. Note that the constituent due to the jth interferer is a function of tunable parameter 

ay- 
Denote by Pj the power of the downlink signal that the base station transmits 
towards the jth remote user. The design objective of strategy can be expressed as follows: 
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Choose the tunable parameters ot y - so that the ratio of the average powers of the 
remote user i and each co-channel interferer j is the same as the ratio of the powers used 
to transmit to the interferer j and the remote user /'. By using the approximation llayll « 



lla/l, 



la/=— 



synth _ j 



p ifi 



„2 _2 



By adding the constraint that a.j must be a positive real number, it follows that, 



2 ? 



when /»-|a,- 1| (7,- > P y - 



10 



«/=- 



2 2 



When P^faf^Pj^j 



° synth _j \ P j\\*j\ 

<jj , then there is no need for null deepening, and therefore a 



value of zero is used for ctj. 



In many systems, the remote user signals have approximately the same average 

power (cr^a^a^^), so that a simplified power based strategy sets the tunable 

parameters Oj as a functions of the ratio of the power used to transmit to the remote user 

to the power for transmission to the remote user, and the ratio of the square of the vector 

2 || II 2 

15 norms of the remote user and the interferer. Mathematically, when P t ||a,- 1| > Pj a j , 



20 



<*/= 



-1 



2 I |2 

Again, when P^ la,- 1| < Pj ay , there is no need for null deepening, which is equivalent 
to a value of zero being used for ay. 

Note that when applied to the covariance modification method, according to this 
strategy 2, the adjustable parameter for any interferer is selected to approximately 
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maintain the same ratio of interferer power to remote user signal power reflected in the 
constituent parts of the modified covariance as the ratio of respective interferer power to 
remote user transmit power used to transmit to the respective interferer and the remote 
user, respectively. 

5 Strategy 3 (applied, for example, to Covariance Modification) 

In addition to the tunable parameter setting method described above using the 
signal injection method as an illustrative example, another feature of the invention is 
determining exactly the parameter setting for an amount of null depth required in null 
deepening when using the any of the inventive null deepening methods. 

10 For sake of simplicity of notation, we'll consider the case where there are two 

dominant co-channel users (e.g., two spatial channels), users i and j. The method can be 
easily extended to more users. The covariance modification method and the strategy 
orthogonalization method will be considered as examples. 

In the linear strategy case, assume a modified linear strategy (spatial or 
15 spatio-temporal) exists for receiving from or transmitting to a desired remote user, say 
the ith, with an interferer j and using parameter oc, to obtain the "null deepened" strategy 
w f - . Similarly, assume that another modified linear strategy exists for receiving from 

or transmitting to an interferer, say the jth user, where j * i, the modification using 
parameter ex,, with w z - and w^-have substantially the same null and gain patterns, except 

20 of course that w,- nulls the interferer and Wj nulls the desired remote user i because when 
communicating with the interferer, the remote user is an interferer. 

According to this third strategy, applicable to the case of the strategy applied in 
the downlink, the parameters ot; and ct^ are chosen to minimize total transmit power, 
subject to the constraints that the signal quality experienced by each co-channel user 
25 (e.g., i and j) in the set including the particular remote user of interest and at least one of 
the one or more interferers present, meets or exceeds some prescribed quality of service. 
Several measures of signal quality may be used, and the signal-to-interference-plus-noise 
ratio (SINR) is preferably used. Using the model of Eq. (13) that 
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Z = a^-f ajsj¥ v, 

where v is the interference (other than user j) and noise, and assuming a modified linear 
strategy w f - and w y - for remote users / and y, respectively, estimates of the SINR of the 

two co-channel users are 



SINR, = 



and 



- H 



SINR, = 



w?a,- 



- H 



Therefore, determining the parameters a £ and aj that minimize the total transmit power 
subject to the SINR being above some threshold SJNR t ^ res h is solving 



argmin {a/ , oc,}{ I + 



w 



} 



10 subject to 

SINR, > SJNR thresh and SINR, > SINR,^, 
where w £ * and w y * are the strategies for user i and j respectively after applying null 
deepening with parameters a t - and cc y - using any of the inventive null deepening method 
embodiments. 

15 As an example, consider the covariance modification embodiment decribed by 

Eq. (11) applied to the weight computation method described by Eq. (7 A). To determine 
the parameter a,, we first apply the well known matrix inversion lemma to the 
combination of Eq. (7 A) and Eq. (1 1) to obtain an expression for w f in terms of the 

quantities , r ZSi = Zs- 4 , a,, a y , a,-, and ay. In particular, 

I ay R zz a 7 a ? R zz r zs7 
20 w,-- k zz r zs/ - 
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Note that except for a,-, and ay, all of the quantities in this expression and in the 
constraint expressions for SINR, and SINRy can either be readily estimated (a,* 2 , aj 1 a /? 
and a,-.) or else are already computed (R~ z , r zs/ , and r^). Hence, the SINR, and SINR, 
constraint equations can be written as 

ii , q,,(^ 2 ) / («i)^0 (14A) 

l=Om=0 

where the constants 'c/ m are defined in terms of a,- 2 , Oy 2 a,-, ay, R zz , r zs/ , r zs p and 
SXNR thresh . 

2 II II 2 

Similarly, the total power cost function P = { || w,- 1| + w j | } can be expressed 

as 

p d iA {a}) 2 +d ia ccj +d i , 3 [ djjaj) 2 +d jt2 aj+dj !3 
d iA («? ) 2 + d i<5 af + 1 d jA (a) ) 2 + d i5 ccj + 1 

Eqs. (14A) and (14B) constitute a constrained nonlinear optimization problem in and 

aj , and can be solved numerically using any of many known nonlinear constrained 

optimization techniques, as would be clear to one of ordinary skill in the art. In practice, 
one interferer tends to dominate among co-channel users, namely the user farthest away. 
Applying null-deepening to this dominant interferer alone accrues most of the gains. In 
this case, the optimization problem is considerably simplified, since there is then only a 

single variable, aj over which the solution needs to be optimized. The constraints are 
then quadratic in a j . 

Note the above description of parameter setting methods is not meant to be 
limiting. It would be clear to those in the art that many other alternate parameter setting 
strategies are possible for the invention. 
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Improved Signature Estimation 

The interferer signatures to use in any of the above embodiments may be known 
signatures, for example, from known records. If the signatures are not known, the 
signatures maybe estimated using any method, including those now described. 

Consider two subscriber units, denoted SU ( - and SUp respectively. When 
communicating with SU/, for example, SUj is an interferer. The null deepening methods 
described above all require known or estimated signatures of the known interferers. An 
improved method for estimating signatures applicable to null deepening is now 
described. 

When the received signals are sampled, say with N samples, the base station 
computes the weight for the ith remote user. Many methods are possible for the weight 
computation, including by the methods described by Eqs. (7 A) and (7B). That is, the 
weight to use for the ith remote may be determined by 



2 

W; = argmin w 



|w H Z- S/ | = (R zz y*r zsi = (ZZH)-i ZsJ 1 , 

or in an improved version, by 

Wf = (ZZH + yl)" 1 ZsJ* , 
where is an estimate, for example the reference signal for the ith transmitted signal, 
and /=1, 2, . . ., Ate, the number of co-channel users. Note that the quantity r zsi = Zs[* is 
computed at the base station for each remote user i in order to form weights to determine 
each of the signal estimates Sj. This observation motivates one known method of 
estimating the signature of the ith remote user, which under the assumption of white 
noise, gives the maximum-likelihood (MLE) estimate when there is only a single 
subscriber unit. That estimate is 



a, =Zs/ y (s i sf f r 1 , (15) 

77 1 

where s { is the reference signal for the ith remote user. The (s,s f - ) term is a scalar 
factor used for normalizing the signature. 
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When other remote users share the same channel, the spatial signature estimate of 
Eq. (15) is contaminated by cross correlations between signatures. This may be seen by 
combining the model of Eq. (4) with Eq. (15), to obtain 



(16) 



where matrix V is 

V = [v(l)v(2)...v(A0]. 

The second term of Eq. (16) is implicitly assumed to be zero in the prior-art signature 
estimation method described by Eq. (15), but this assumption is inaccurate for short data 
sequences as encountered in practical systems, such as in the wireless local loop (WLL) 
system base station and the mobile PHS base station in which the method of present 
invention is implemented. 

Another method for signature estimation takes these cross-correlation terms into 
account/Let S is the matrix with the Ns reference signals as the rows. That is, 



S = 



si 



Let A be a matrix having the Ns signature values as columns. That is, A = [a t a 2 . . . a^J. 
Then the MLE solution is 

A = arg min ||Z - AS|| 2 =ZS H (SS H )" 1 (17) 
A 

where A is the matrix having the signatures estimates as columns, 

A = [a 1 a 2 ■■• &n s ]. (18) 
See for example, A. Lee Swindlehurst, S. Daas, J. Yang: "Analysis of a Decision 
Directed Beamformer, IEEE Trans, on Sig. Proc. Vol. 43, No. 12, pp. 29202927, 
December, 1995. 



For Ns=2 remote users, the estimates of Eq. (17) can be expressed as 
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z sPhll 2 -^(s 2 sh 

hi 2 w 2 -«rf 2 



(19) 



and 



a 2 =Z 



(20) 



respectively. 

While the estimator of Eq. (17) gives an accurate signature estimate, it is 
computationally complex. The third signature estimation method described herein is an 
iterative method which is computationally simpler than that of Eq. (17) but provides for a 
more accurate estimate than that of Eq. (15). 

The method is first described for the case that the base station is (or is assumed to 
be) handling two spatial channels {i.e., two remote users communicating with the base 
station on the same carrier and in the same time slot), or that there is one spatial channel 
and a known interferer. The signature estimate for the second remote user (the interferer) 
can then be computed as follows: 

Step 1: Compute the an estimate of aj using any known method, for example, that 



of Eq. (15): 



a 1 =Zs 1 H (s 1 s 1 H r 1 
where is the reference signal for SU1. 



(21) 



Step 2: Subtract from the received signal an estimate of the contribution of the 



first remote user: 



Z = Z-a,s 



(22) 



Step 3: Compute the improved estimate of a 2 using the received signal with the 



first user estimated contribution removed: 
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a 2 =Zs?(s 2 s?)- 1 



(23) 



where s 2 is a reference signal obtained for the second remote user. That is, 
one needs to have a reference signal for the interferer towards which one is 
directing the null. 

Since the correction term subtracted from Z in Eq. (22) is itself an approximation, 
iteratively repeating Eqs. (22) and (23) may improve the estimate further. For the case in 
which Eq. (15) is used in step 1 above, substituting Eqs. (21) and (22) into Eq. (23), one 
obtains 



Eqs. (24) and (20) differ only by a constant, so the recursive method in this case 
converges to the MLE solution (within a scale factor) in one iteration. 

This recursive method can be extended to incorporate an arbitrary number of 
remote users. A method for Ns=3 is as follows. To estimate the signature for some 
additional remote user, denoted SU3, with three remote users Ns=3: 

Step 1: Compute a 2 using any known method, for example the method of Eq. (21) 
(applied however to SU2). 

Step 2: Form the decorrelated estimate which removes the SU2 contribution 



a 2 =Z^-s 1 H (s 1 sP)- 1 s 1 }?(s 2 s^r 1 




sF( Sl s??) 



(24) 



a 2 =Zs?(s 2 s?r 1 



(25), 



where Z is given in Eq. (22), but applied for SU2. 



Step 3: Iterate the estimate of a t with the value of a 2 just found, 



5 1 1 =(Z-a 2 s 2 )s 1 H (s 1 sP)- 1 



(26) 



Step 4: Compute the estimate of the desired SU 3 signature, 



a 3 = (Z - a/s! - a 2 s 2 )sf (s 3 sf ) 



(27) 



where S3 is a reference signal for the user SU3. 



• 



40 

Estimation of signatures for two remote users, SU1 and SU2, with Ns~3 results 
from appropriately changing the indices in this procedure. 

A generalization of the iterative method is as follows for the general case of Ns 
signatures to be estimated: 

Step 1: Assume some initial estimates. For example, assume the starting values 
are those given by Eq. (15), i.e., a, = Zs^ 1 (s/sj 1 ) _1 for i = 1, . . Ns. 

Step 2: Obtain an estimate for each signature, while keeping the values of the 
remaining Ns-l signatures to be the most recently determined estimate 
values, each of these fixed most recently determined estimate values initially 
being the initial value from step 1, and later being either the values 
determined when this Step 2 was previously applies to another signature in 
the same iteration, or, for those signatures for which no estimate has yet been 
obtained during the current iteration, the estimate from the most recent 
previous iteration. Denoting the jth user's estimate by a,, this jth signature 
estimate is obtained as the solution to 

&j = arg min ||Z - AS|| 2 , 

while the other Ns-l values are so fixed, where again A = [aj a 2 . . . a^]. 
Repeat this procedure for ally = 1, 2, . . Ns, replacing the value at which any 
signature estimate is fixed to be the most recently obtained estimate value. 

Step 3: Repeat Step 2 iteratively until convergence is reached. 

To further improve signature estimation, in another aspect of this invention, each 
desired signature estimate &j is singular value decomposition-averaged (S VD-averaged) 
over successive signal streams {i.e., bursts) to reduce the contribution of the noise-related 
error term in Eq. (15). One method of performing SVD averaging comprises stacking up 
the signatures to be combined as the columns of a matrix, then performing the SVD of 
the resulting matrix. The SVD averaged value is then taken as the left singular vector 
corresponding to the largest principal value. An alternate SVD-averaging embodiment 
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uses a weighted moving window (last few signatures) SVD average which favors recent 
data. Other combining methods also may be used. 

To see how having a better signature estimate improves the null deepening, 
consider the case of two users, user SU, and user SUy, with SU t being the desired user 
and SUy the interferes and let 

ay = ay + |i a,- + e 

be the estimate of the interference signature ay, where e is the estimation error orthogonal 
to a,-, the signature of user SU ( , and (|l a,) is the estimation error in the direction of a,-. 

The null deepening methods described herein direct a null in the direction of a,- 
and a gain in the direction of a,-. Thus, if we want a good null in the direction of ay, we 
must have both \i very small, and a good null in the direction of 8 (even if e is 
comparatively large). One case for which the term |i may be undesirably large occurs 
when the desired remote user signal is present in the received data and when the 
signature a,- is estimated using Eq. (15), i.e., using cross-correlation with the reference 
signal sj. A smaller value of allows a deep null to be directed in the direction of a,-. The 
term \i can be greatly reduced by estimating the interferer signature a,- using Eq. (17) for 
joint maximum likelihood estimation on the reference signals corresponding to both the 
remote user and the interferer, or by using the iterative method. In each case, a deep null 
may then be directed in the direction of a,. 

Results of Simulations 

To see how having a better signature estimate improves the null deepening, 
consider the case of two users, user SU, and user SUp with SU, being the desired user 
and SUj the interferer, and let 

be the estimate of the interference signature ay, where e is the estimation error. 

Assume also that for the purpose of this performance analysis, the received signal 
may be modeled as 
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Z= a; s/ + ay sy + v, 
where v represents temporally white noise and interference. 

For the signature estimation method described in Eq. (15) and for the improved 
method described in Eq. (17), the signature error covariance E[e e H ] is, up to a constant 
5 of proportionality, approximately equal to E[v v H ] + a £ a,-* 1 , where for the method of 
Eq. (15), ^ is proportional to the SNR of SU/ and for the method of Eq. (17), |l is nearly 
zero. 

In one embodiment of the invention, a weight w f . is computed that is close to the 
MMSE weight w,- such that w . is nearly orthogonal to ay. Thus, 

10 0 = w ■* ay = w • ay + \v, H e. 

If we choose to measure null depth by the expected squared dot product 

ATD = E[lwf t ayl 2 ], 
then the null depth approaches zero (the depth of a perfect null) if 

E[l w J 1 eP] = E[ w} 1 E[e e H ] w . ] 

15 = E[ wf 1 E[v vH] w . ] + ME[ w J 1 a/ a/ H w . ] 

also approaches zero. The term E[ w - 1 E[v v H ] w f . ] in the right hand side of the above 

equation tends to be small for any interference-mitigating weight \v i (e.g., for a weight 
determined using a reference signal of the desired user signal, such as the methods of 
Eq. (7 A) or Eq. (7B)). The |iE[ w J 1 a/ a/ H w f ] term, however, is small only if |i is small, 

20 because the weight w f - delivers power towards the signature a/ of the desired user SU;. 

Thus the method of Eq. (17) for signature estimation — which has |X close to zero — will 
often yield improved null depth compared to the method of Eq. (15). 

Simulations were used to evaluate the performance of several embodiments of the 
present invention. 
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The uplink signal was formed by adding three equal-power simulated subscriber 
unit signals, each consisting of 50 symbol-spaced samples of a randomly generated tc/4 
DQPSK baseband waveform. The simulation used a 12-element antenna array for both 
transmission and reception. Uplink spatial signatures were drawn randomly and 
uniformly from the complex 12-sphere; uplink and downlink signatures were assumed 
equal. Additive white Gaussian noise was added to the uplink signal to set the received 
SNR equal to 12 dB, where SNR is measured at the output of an ideal interference-free 
linear beamformer. 

Performance was computed by Monte-Carlo averaging over 1000 independent 

runs. 

Signature estimation was performed using both the single-user maximum- 
likelihood method of Eq. (15) and the two-user maximum-likelihood method of Eq. (17). 
Signatures were averaged using SVD averaging over Navgs = 1, 2, 4, 8, and 16 bursts. 
The accuracy of signature estimation was evaluated using the dot product 



c = 



a 5S 



mm 



Three nulling strategies were computed: the baseline MMSE strategy of Eq. (7 A), 
a null-deepened strategy computed using signal injection method with oc, = 1000, and a 
null deepening strategy using weight orthogonalization. Null depth was measured by the 
expected dot product 

ND = ±- 



The simulated received data used to compute the signature estimate was reused to 
compute the null-deepened strategy. 

The null depth WD MMSE of the baseline prior-art MMSE method was compared 
to the null depths M> s ig_i n j and #A;ig_inj °f the inventive signal injection and inventive 
orthogonalization methods, respectively, using the quantities 
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rfB sig_inj = 20 log 10 (Atf>MMSE / ^sig_inj). 

and 

^orthog = 20 ^glO^MMSE 1 ^orthog). 
respectively. Median results over 1000 runs are tabulated below. 

5 



Navgs 


c (eq 15) 


c (eq 17) 


*^sig_inj 


^orthog 


1 


0.9792 


0.9869 


0.86 dB 


0.78 dB 


2 


0.9896 


0.9930 


3.80 dB 


3.94 dB 


4 


0.9942 


0.9965 


6.76 dB 


6.75 dB 


8 


0.9971 


0.9984 


9.94 dB 


9.93 dB 


16 


0.9987 


0.9992 


12.5 dB 


12.6 dB 



As can be seen from the table, the results were almost identical for the two 
embodiments of the invention. Also, as can be seen from the table and from Fig* 16, 
doubling the number of bursts over which the signature estimate is averaged increases the 
null depth improvement by roughly 3 dB. 

10 The simulations also confirmed that the depth of the null directed towards the 

second interferer (i.e., the interferer towards which no deepening was performed and no 
signature estimation carried out) was not significantly affected by the signal injection 
method or orthogonalization embodiments of the inventive method. 

Other Embodiments 

15 In another embodiment, an explicit signature estimate is not needed, and instead 

the uplink weights for two co-spatial remote users are combined to construct weights 
with improved nulling performance, while again substantially preserving the rest of the 
null and gain pattern. 

Note that the apparatus and method of the present invention does not need to 
20 explicitly identify any new interferer waveforms. Note also that the various embodiments 
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of the invention can be used for smart antenna processing on the downlink, or for smart 
antenna processing on the uplink. That is, the embodiments can be used on the downlink 
for smart antenna processing signals that are received at the antenna elements of an 
antenna array which is part of a wireless receiver, or on the uplink for smart antenna 
processing a signal to be transmitted to form signals for transmission by the antenna 
elements of an antenna array which is part of a wireless transmitter. In general, the term 
"wireless station" therefore may be used for: (i) a receiver having an array of antenna 
elements; (ii) a transmitter having an array of antenna elements; or (iii) a transceiver 
having an array of antenna elements for both reception and transmission. 

Several aspects of the invention described herein were described implemented as 
programs run on one or more DSP devices. Given sufficient economic incentive, DSP 
functionality, including DSP programs, may be incorporated into special purpose 
hardware, for example as part of an application specific integrated circuit (ASIC) or as 
part of a very large scale integrated circuit (VLSI). DSP functionality may also be met by 
other processors, for example a general purpose microprocessor. In addition, a DSP 
device running a program may be converted into a special purpose piece of hardware. 
Thus, the terms processor, digital signal processor, DSP, and DSP device as used herein 
include these equivalent alternatives. 

The terms and expressions which have been used herein are used as terms of 
description and not of limitation, and. there is no intention in the use of such terms and 
expressions of excluding equivalents of the features shown and described, or portions 
thereof, it being recognized that various modifications are possible within the scope of the 
invention claimed. 



